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Abstract 

The Skyrme-Faddeev model is a modified sigma model in three-dimensional space, 
which has string-like topological solitons classified by the integer-valued Hopf charge. 
Numerical simulations are performed to compute soliton solutions for Hopf charges up 
to sixteen, with initial conditions provided by families of rational maps from the three- 
sphere into the complex projective line. A large number of new solutions are presented, 
including a variety of torus knots for a range of Hopf charges. Often these knots are 
only local energy minima, with the global minimum being a linked solution, but for 
some values of the Hopf charge they are good candidates for the global minimum 
energy solution. The computed energies are in agreement with Ward's conjectured 
energy bound. 
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1 Introduction 



Over thirty years ago Faddeev suggested [I] that in three-dimensional space the 0(3) 
sigma model, modified by the addition of a Skyrme term, should have interesting string-like 
topological solitons stabilized by the integer-valued Hopf charge. Ten years ago substantial 
interest was generated by the first attempts at a numerical construction of such solitons 
[5J E] and the suggestion that minimal energy solitons might take the form of knots [5] . This 
conjecture has been confirmed by numerical results, which demonstrate that the minimal 
energy soliton with Hopf charge seven is a trefoil knot p]. Substantial numerical investiga- 
tions by several authors [H El [H El EE] has produced a comprehensive analysis of solitons 
with Hopf charges from one to seven, and it appears that the global energy minima have 
now been identified for these charges, together with several other stable soliton solutions 
that correspond to local energy minima. 

For Hopf charges five and six the minimal energy solitons form links, but so far the 
charge seven solution remains the only known knot, even including local energy minima. It 
is therefore currently unknown whether this charge seven trefoil knot is the only knotted 
solution or if there are many other knots, perhaps of various types, which arise for higher 
Hopf charges. It is this question which is addressed in the present paper. Solitons with Hopf 
charges up to sixteen are constructed numerically and it is found that a variety of torus 
knots exist at various Hopf charges. Often these knots are only local energy minima, with 
the global minimum being a link, but for some values of the Hopf charge it appears that a 
knot is a good candidate for the global minimum energy solution. 

One of the difficulties in extending previous numerical studies to higher Hopf charges is 
that, even at low charges, there are a number of local energy minima with large capture 
basins. This makes it difficult to fully explore the landscape of local energy minima and 
hence determine the global minimum, with the severity of the problem generally increasing 
as the charge increases. To overcome this problem an analytic ansatz is employed which uses 
rational maps from the three-sphere into the complex projective line. Although this ansatz 
does not provide any exact solutions, it does allow the construction of reasonable initial 
conditions for any torus knot for a large range of charges, plus a wide selection of links and 
unknots. It is therefore possible to start at a variety of locations in field configuration space, 
in particular where one suspects that a local energy minimum might be close by, and explore 
the energy landscape around this point. Not only is this a fruitful approach for finding new 
solutions, but it can also be used to provide strong evidence that a solution of a particular 
type does not exist, by starting with an initial condition approximating this configuration 
and demonstrating that it changes dramatically under energy relaxation. 

An analysis is made of the new knotted and linked solutions and some understanding is 
obtained regarding the charges at which particular knots are likely to exist. 
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2 The Skyrme-Faddeev model and low charge solitons 



The Skyrme-Faddeev model involves a map : M 3 i— > S 2 , which is realized as a real three- 
component vector cf) = (0x, 02, 03), of unit length, <f> ■ (ft = 1. As this paper is concerned only 
with static solutions then the model can be defined by its energy 

E = / d l( f> ■ d i( j> + ~(c^ x d ]( t>) ■ (<9,0 x d 3 <t>) d 3 x, (2.1) 

where the normalization is chosen for later convenience. The first term in the energy is 
that of the usual 0(3) sigma model and the second is a Skyrme term, required to provide a 
balance under scaling and hence allow solitons with a finite non-zero size. 

Finite energy boundary conditions require that the field tends to a constant value at 
spatial infinity, which is chosen to be 0(oo) = (0, 0, 1) = e 3 . This boundary condition com- 
pactifies space to S 3 , so that the field becomes a map (f) : S 3 i— > S 2 . Such maps are classified 
by n 3 (S 2 ) = Z, so there is an integer- valued topological charge Q, the Hopf charge, which 
gives the soliton number. Unlike most theories with topological solitons, the topological 
charge is not a winding number or degree of a mapping. Rather, it has a geometrical inter- 
pretation as a linking number of field lines, as follows. Generically, the preimage of a point 
on the target two-sphere is a closed loop (or a collection of closed loops) in the three-sphere 
domain obtained from the compactification of M 3 . Two loops obtained as the preimages of 
any two distinct points on the target two-sphere are linked exactly Q times, where Q is the 
Hopf charge. 

The Hopf charge can be written as the integral over the three-sphere domain of a charge 
density, but this density is non-local in the field </>, as the following construction demon- 
strates. Let uj denote the area two-form on the target two-sphere and let F = 4>*u be its 
pull-back under (f> to the domain three-sphere. The triviality of the second cohomology 
group of the three-sphere implies that F is an exact two-form, say F = dA. The Hopf charge 
is then given by integrating the Chern-Simons three-form over the three-sphere as 



The energy bound 



/ 3 x 3/8 

E>cQ 3/ \ where c = f — J ^ 0.534 (2.3) 

has been proved [T2l 19]. and it is known that the fractional power is optimal [TO], though it is 
expected that the above value of the constant c is not. Motivated by a study of the Skyrme- 
Faddeev model on a three-sphere with a finite radius, together with an analogy with the 
Skyrme model, Ward has conjectured [13] that the above energy bound holds with the value 
c = 1 (hence the choice of normalization factor for the energy in (12.11) ) , but this has not been 
proven. As discussed shortly, the current energy values known for low charge solitons are in 
good agreement with Ward's conjectured energy bound, and later it will be demonstrated 
that the energies for higher charge solitons are too. 
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The simplest examples of Hopf solitons are axially symmetric and their qualitative fea- 
tures may be described as follows. First, consider the model ( 12.1 j) defined in two-dimensional 
space, so that the boundary conditions now result in a compact ification of space from M 2 to 
S 2 . The field is then a map <fi : S 2 i— > S 2 and is classified by an integer m £ Z = ^(S* 2 ), which 
is the usual winding number familiar from the 0(3) sigma model itself. Strictly speaking, 
in order to have two-dimensional solitons with a finite size in this model then a potential 
term must also be included, upon which the two-dimensional solitons are known as baby 
Skyrmions [UJ, but this aspect is not crucial for the discussion here. Returning to the 
three-dimensional model, then a toroidal field configuration can be formed by embedding 
the two-dimensional soliton in the normal slice to a circle in space. The two-dimensional 
soliton has an internal phase and this can be rotated through an angle 27m, as it travels 
around the circle once, where n £ Z counts the number of twists. A field configuration of 
this type has Hopf charge Q = nm, and it will be denoted by A n , m , so the first subscript 
labels the number of twists and the second is the winding number of the two-dimensional 
soliton forming the loop. 

The minimal energy Q — 1 soliton is of the type Ai,i and has been studied numerically by 
a number of authors 0, El E El EH] using different algorithms. Taken together, these results 
suggest that the charge one soliton has an energy of around E = 1.21, and this is part of 
a general pattern where the minimal energy soliton exceeds Ward's conjectured bound by 
around 20%. The numerical computations reported in this paper produce a charge one energy 
of E = 1.204, and hence appear to underestimate the energy by an amount of the order of 1%. 
For low charge solitons this accuracy could be improved by increasing the resolution of the 
numerical grid, but it is computationally too expensive to increase the resolution for all the 
large number of simulations reported here for large charges. Moreover, although the energy 
is a slight underestimate for a given soliton solution, the relative energies between two given 
solutions are much more accurate than this (changing grid resolutions and sizes suggests an 
accuracy of around 0.1%) and it is more important to preserve the relative energies than the 
absolute energy of any single solution. 

The details of the numerical computations are similar to those of Ref. [I], and full details 
of the related numerical code used to study Skyrmions can be found in Ref. [2]. Briefly, a finite 
difference scheme is employed with a lattice spacing Ax = 0.1 on a grid containing (151) 3 
lattice points, with the field fixed to the vacuum value <fi = e 3 on the boundary of the grid; 
this is preferable to other methods of dealing with the finite volume simulation domain, such 
as the one used in Ref. [I] which led to a more substantial underestimate of energies, though 
again relative energies had a high accuracy. The energy minimization algorithm proceeds by 
evolving second order in time dynamics and periodically removing kinetic energy from the 
system whenever the potential energy of the system begins to increase. It is based on the 
approach employed in Ref. [I], but is computationally more efficient since the second order 
dynamics is derived from the kinetic term of the sigma model only, rather than the full 
kinetic term determined from the Lorentz invariant Lagrangian associated with the energy 
(12. ip . This subtle difference leads to a more efficient algorithm since the leading term in the 
evolution equations is now diagonal, so a costly matrix inversion can be avoided. 
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A sensible definition of the position of a soliton is to identify where the field is as far 
as possible from the vacuum field. A Hopf soliton is therefore string-like, since the soliton's 
position is defined to be the closed loop (or collection of loops) corresponding to the preimage 
of the point <fi = —e% = (0,0, —1), which is antipodal to the vacuum value on the target 
two-sphere. To visualize a field configuration this position curve is plotted, though for clarity 
a tube around this position is displayed, given by an isosurface of the form 3 = — 1 + v, 
where generally the value v = 0.2 is chosen. If a pictorial representation of the Hopf charge 
is also required then a second preimage curve can also be plotted (again as a tube), and 
the linking number of these two curves inspected. There is no natural choice for this second 
preimage value (e 3 is not a very useful choice) and in this paper when a second preimage 
curve is displayed it corresponds to the point 4> = (y/2fi — /i 2 , 0, —1 + jj) where /i = 0.1. In 
what follows this curve will be referred to as the linking curve, and since \i is relatively small 
then the linking curve remains reasonably close to the position curve, making it fairly easy 
to inspect the linking number. 

In Figure [I] the position (light tube) and linking (dark tube) curves are displayed for 
the known lowest energy solitons with Hopf charges from one to seven. Given the large 
number of numerical simulations performed by different groups it seems reasonably certain 
that these are the minimal energy solitons for these charges. The energies of these solutions 
will be discussed later, and as mentioned above they are generally around 20% above Ward's 
conjectured bound, but for now it is the structure of the solutions that is of primary interest, 
and this is reviewed below. 

The solutions with charges one and two are both axially symmetric, and are of the type 
A.\,\ and .4.2,1 respectively [51 E], with one and two twists of the linking curve around the 
position curve clearly visible. There is a solution of the form A\p. p3J [HI [7] but it has 
an energy around 13% above that of the minimal charge two soliton. This solution has 
the interpretation of two Q = 1 solitons stacked one above the other, preserving the axial 
symmetry. 

The charge three soliton is basically of the type A3,i but the axial symmetry is broken as 
the position curve is bent [1], so it will be denoted by .4.3,1 to emphasize that the structure 
is deformed. There is an axial solution of the type *4 n ,i for any integer n, but for n > 2 this 
solution is unstable to a coiling instability [1]; a well-known phenomenon in other settings 
such as twisted elastic rods. The Q = 4 soliton is of the type ^2,2, [HI EES], and may be 
thought of as two Q = 2 solitons stacked one above the other. For both Q = 4 and Q = 5 
there are also bent solutions PQ, similar to the charge three soliton, and therefore denoted 
.4.4,1 and .45,1, but unlike the Q = 3 case, both these solutions are only local energy minima, 
with energies around 2% and 4% above the global minimum, respectively. The minimal 
energy charge five soliton provides the first example of a link, that is, a solution in which 
the position curve contains two or more disconnected components. It consists of a charge 
one soliton and a charge two soliton which are linked once [8], with the result that each of 
these solitons gains an additional unit of linking number to make a total of five. To denote a 
link of this type the notation £\'^ will be used, where the subscripts label the charges of the 
components of the links and the superscript above each subscript counts the extra linking 
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Figure 1: The position (light tube) and linking (dark tube) curves for the known lowest 
energy solitons with Hopf charges 1 < Q < 7. 



number of that component, due to its linking with the others. Hence the total charge is the 
sum of the subscripts plus superscripts. The charge six soliton is similar to that of charge 
five, but now both components of the link have charge two pQ, so using the above notation 
it is written as £^2- 

Finally, the charge seven soliton is the first (and so far only) example of a knot pp. The 
position curve is a trefoil knot, which has a self-linking or crossing number of three (the 
crossing number of a knot is not an invariant and its use in this paper refers to the minimal 
crossing number over all presentations). An examination of the linking curve in Figured] 
confirms that it twists around the position curve four times as the knot is traversed, so the 
Hopf charge is indeed seven, being the sum of the crossing number plus the number of twists. 
Any field configuration of a trefoil knot will be denoted by /C3 2, which refers to the fact that 
the trefoil knot is also the (3, 2)-torus knot (see the following section for more information 
regarding torus knots). Note that this notation does not display the Hopf charge of the 
configuration, which is determined by the twist number as well as the crossing number, but 
this should not cause any confusion in what follows. Of course, for knots and links, as well 
as for unknot configurations, the precise conformation of the position curve is important in 
determining the energy, not just its topological type, but its knotted or linked structure is 
an important feature in classifying the solution. Thus, for example, the Q = 7 trefoil knot 
displayed in Figure [T] is not as symmetric as its typical knot theory presentation, and it is 
significant that the energy is lowered by breaking the possible cyclic C3 symmetry. 

Having reviewed the known results for low charge solitons and introduced the notation 
used to label their structural type, it is time to turn to the main questions addressed in 
this paper. Given the variety of solutions which appear, even at low charges, it is difficult 
to make a confident prediction of the kind of behaviour that might arise at higher charges. 
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In particular, an interesting open question is whether more knot solitons appear at higher 
charge, and if so what types of knots arise and are they local or global minima? It could 
have been the case that the charge seven trefoil knot was the only knot soliton, though in 
this paper it will be shown that this is far from true. Also, it would be useful to have at least 
a qualitative understanding of the fact that Q = 7 is the lowest charge at which a knot first 
appears, and desirable to be able to estimate the charges at which other knots might exist. 

In this section it has been mentioned briefly that at each charge there are often local 
minima in addition to the global minimum, and it is expected that the number of local 
minima generally increases with the charge. As noted in the introduction, these local minima 
can have large capture basins, making it difficult to fully explore the space of solutions. For 
any relaxation computation this makes the choice of initial conditions a crucial issue. In 
particular, it is not obvious how to construct a variety of reasonably low-energy initial 
field configurations for all Hopf charges. One family of configurations are the axial fields 
A n ,i, which are all unstable for n > 2. However, even for charges as low as four, initial 
conditions which use small perturbations of these unstable solutions get trapped in local 
minima. One successful approach for low charges involves constructing links by hand [3 [8], 
using a numerical cut-and-paste technique where various numerical field configurations are 
sewn together in different parts of the simulation grid to create a linked field with a given 
Hopf charge. Although this approach was the first method to successfully yield the global 
minima at charges four and five, it becomes more cumbersome for larger charges and is not 
applicable for creating knot initial conditions; it is also not very elegant mathematically, but 
this is perhaps not a serious criticism. Another set of initial conditions that have been used 
[8] are based on the fields A n , m , with m > 1, though these fields tend to have relatively 
large energies and therefore require quite long simulation times (except for the low charge 
examples of Ai } 2 and ^2,2) • In the following section an analytic ansatz is presented that yields 
reasonable field configurations for a variety of charges and describes all torus knots plus a 
wide range of links, in addition to fields of the type A ntm . The key ingredient of the ansatz 
involves a rational map from the three-sphere into the complex projective line. The ansatz 
allows the construction of initial conditions corresponding to a variety of locations in field 
configuration space; in particular, configurations can be created which are reasonably close 
to suspected local energy minima. This approach substantially reduces the computational 
effort required to explore the energy landscape and makes it feasible to study knots and links 
for quite large charges. 

3 Rational maps, torus knots and links 

Torus knots are classified by a pair (a, b) of coprime positive integers with a > b. A knot 
is a torus knot if it lies on the surface of a torus, in which case the integers a and b count 
the number of times that the knot winds around the two cycles of the torus. The simplest 
example is the trefoil knot, which is the (3, 2)-torus knot and has three crossings. In the 
standard knot catalogue notation it is 3i, where the number refers to the crossing number 
of the knot and the subscript labels its position in the knot catalogue. Other torus knots 
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which will be of interest in this paper include the (5, 2)-torus knot, also known as Solomon's 
seal knot, which has five crossings (knot 5i in the catalogue), and the (4, 3)-torus knot which 
has crossing number eight (knot 819 in the catalogue). In general the (a, b) -torus knot has 
crossing number C = a(b — 1). 

Denote by K a ^ any field configuration in which the position curve is an (a, 6)-torus knot. 
Such a field configuration will have Hopf charge Q = C + T, where C = a(b — 1) is the 
knot crossing number and the integer T counts the number of times that the linking curve 
twists around the position curve (only situations where the orientation is such that T is 
non- negative will be of relevance here). The main aim of this section is to construct fields 
of type /C a> 6 for a large range of charges Q. 

The first step is to recall the standard description [3J of a torus knot as the intersection 
of a complex algebraic curve with the three-sphere. Consider (Zi, Zq) 6 C 2 and the unit 
three-sphere S 3 C C 2 given by \Zi\ 2 + \Z \ 2 = 1. Then the intersection of this three-sphere 
with the complex algebraic curve Zf + Zq = is indeed the (a, 6)-torus knot. 

To use the above description to produce a field configuration, the spatial coordinates 
(xx, a?2, £3) G M 3 are mapped to the unit three-sphere via a degree one spherically equivariant 
map. Explicitly, 

sin f sin J 

(Z 1 ,Z ) = ((xi + ix 2 )—^—,cosf + i—^—x 3 ), (3.1) 

where r 2 = x\ + x\ + x\ and the profile function f(r) is a monotonically decreasing function 
of the radius r, with boundary conditions /(0) = tt and /(oo) = 0. The precise form of this 
function will be specified shortly. 

A Riemann sphere coordinate W, is used on the target two-sphere of the field 0, deter- 
mined via stereographic projection as 

W=^±^. (3.2) 

1 + 03 

With these definitions, then constructing a map <fi : M 3 1— > S 2 is equivalent to specifying 
W(Zi, Z ), that is, W : S 3 1— > CP 1 , a map from the three-sphere to the complex projective 
line. This map will be taken to be a rational map, which means that W = p/q, where both 
p and q are polynomials in the variables Z\ and Zq. 

For an (a, 6)-torus knot the mapping is chosen to have the form 



„ 17 a 7/3 



where a is a positive integer and (3 is a non-negative integer. A map of this form has 
the correct boundary conditions at spatial infinity, since as r — > 00 then (Z 1; Z ) — * (0,1) 
and, because a > 0, this gives W — * 0, which is the stereographic projection of the point 
(ft = (0, 0, 1) = 63. Furthermore, the position curve is the preimage of the point — e3, which 
corresponds to W = 00, and hence is given by q = = Z\ + Z\. As the coordinates (Z±, Zq) 
automatically lie on the unit three-sphere then the position curve is the (a, 6)-torus knot, for 
all values of a and (5. 
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The Hopf charge of the map determined by (13. 3p is given by the cross ratio 



Q = ab + (3a. 



(3.4) 



To see this note that a rational map 



W : S 3 



i— ► 



c: 



:p- 



given by W{Z U Z Q ) 



(3.5) 



has a natural extension to a map 



(p, (/):B 4 CC 2 h 



C 2 , given by (p(Z u Z ), q(Z 1 , Z )). 



(3.6) 



Here B 4 denotes the 4-ball with boundary S 3 , and is obtained by replacing the constraint 
\Zi\ 2 + \Zq\ 2 = 1 by the inequality \Zi\ 2 + \Zq\ 2 < 1. The map (p, q) is between manifolds of 
the same dimension and has a degree, defined by counting preimages of a generic point of 
the target space, weighted by the signs of the Jacobian. As all the other mappings involved 
in the ansatz have degree one, and the standard Hopf map (p,q) i— > p/q is employed, then 
the Hopf charge Q of the combined mapping is equal to the degree of the mapping (p,q). 
Furthermore, as the mapping (p, q) is holomorphic then the sign of the Jacobian of all generic 
preimage points is positive, so the degree is simply the number of preimages. 

For the mapping (p, q) given by (13.31) take the generic point in target space to be (e, 0), so 
the Hopf charge is the number of solutions of the equation (Z^Zq , Z" + Zq) = (e, 0). Writing 
Z in terms of Z\ from the first component of this equation and substituting into the second 
component yields the polynomial equation Zi b+l3a = (— l^e 6 , which clearly has Q = ab + (3a 
solutions, as stated. 

For any (a, 6)-torus knot the above ansatz can be used to construct a field configuration 
of the type K, a ^ with Hopf charge Q, providing there are integers (a, (3), with a positive 
and (3 non-negative, for which Q = ab + (3a. For example, in the case of the trefoil knot 
then (a, b) = (3, 2) hence any Q > 3 can be obtained. In fact a trefoil knot with Q = 3 
can also be obtained by interchanging the roles of Z\ and Z , but a field of this type is not 
relevant, as will be made clear later. Note that for some values of Q, and fixed (a, b), there 
are multiple solutions for (a, (3). These multiple solutions produce qualitatively similar fields, 
in that both describe the same (a, 6)-torus knot position curve, and the linking curve twists 
the same number of times around the position curve; though the distribution of twist may 
vary. In examples where such multiple solutions were used as initial conditions, the energy 
relaxation algorithm produced identical final results for all choices. 

Recall that a profile function f{r) needs to be specified, satisfying the boundary condi- 
tions /(0) = it and /(oo) = 0. In fact, these are the correct boundary conditions in an infinite 
domain but in the finite domain used for numerical simulations the last boundary condition 
needs to be replaced by / = at the boundary of the numerical grid. As the ansatz is 
only used to provide initial conditions then most reasonable monotonic functions will suffice. 
The simulations discussed in this paper were performed on a cubic grid of side-length 2L, 
that is, — L < xi < L, where 2L = 15. A simple linear profile function was used, given by 
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Figure 2: The position (light tube) and linking (dark tube) curves for initial conditions 
created using the rational map ansatz: (a) Q — 7 trefoil knot; (b) Q — 6 link £ 2 ' 2 ■ 



f(r) = tt(L — r)/L, for r < L and zero otherwise. For a given rational map one could aim to 
minimize the energy of the ansatz over all profile functions, but this has not been attempted 
for reasons discussed shortly. 

As an illustration of the use of the ansatz, consider an approximation to the minimal 
energy Q = 7 trefoil knot. To construct a field K^^ with Q = 2a + 3/5 = 7, requires 
(a,/3) = (2, 1), resulting in the map W = ZfZ /(Zf + Zq). The field generated from this 
rational map is displayed in Figure^, where both the position and linking curves are shown. 
It is clear from this figure that the field has similar qualitative features to the minimal energy 
soliton, in that it is a trefoil knot with four twists. The main qualitative difference is that 
the ansatz produces a more symmetric field, having a cyclic C3 symmetry. This symmetry 
is obvious from the form of the rational map, as a rotation by 120° around the a^-axis 
corresponds to the transformation Z\ 1— > e 2m / 3 Zi, upon which the rational map changes by 
only a phase; which is simply an action of the global 0(2) symmetry on target space. 

Using the above field as the initial condition for an energy relaxation simulation yields the 
minimal energy soliton very quickly, and certainly requires far less computational resources 
than computing the solution from a perturbed axial configuration. Note that, in practice, 
the cyclic C3 symmetry is slightly broken by the cubic numerical grid, and in particular 
its boundary, so an explicit symmetry breaking perturbation is not required. In examples 
discussed later, where initial conditions have a cyclic C4 symmetry, then the configuration 
is created with a slight displacement from the centre of the grid to again slightly break any 
exact symmetry. 

As mentioned above, the profile function could be optimized to find the minimal energy 
configuration, for a given rational map, but this has not been attempted. As with the 
above example, generally the ansatz is more symmetric than the relaxed solution, so it is 
not expected that a good estimate of the energy can be obtained in this way. It might be 
possible to improve the approximation to a quantitative level by minimizing over families 
of rational maps, which include symmetry breaking terms, rather than using the rational 
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maps in this paper; chosen to be the simplest with the correct qualitative features. In this 
case it would be worthwhile to also minimize the profile function, but as the minimization 
problem couples the rational map and the profile function together then this is not a simple 
numerical task. In fact, this numerical problem appears to have about the same level of 
difficulty as performing energy relaxation via full field simulations. As the energy relaxation 
code produces a solution very quickly and efficiently from the simple ansatz, then there 
appears little need to try and improve upon it. The profile function could be adjusted to a 
minimal extent, for example, by including parameters to set the scale and thickness of the 
knot, but again this has not been necessary. 

So far only knot initial conditions have been considered, but the ansatz can also be used 
to construct a variety of unknots and links by choosing suitable rational maps. For example, 
axial fields of the form A n , m are generated by a rational map W = Z™/Z™. 

Links correspond to rational maps in which the curve determined by the denominator 
is reducible. The torus knot maps with (a, b) coprime degenerate to links if (a, b) are not 
coprime. As an example, links of the type £^ (which have Hopf charge Q = 2n + 2) are 
derived from the rational map 

Zl-Zl 2(Z 1 -Z ) + 2(z 1 + z y 

where the partial fraction decomposition reveals the charges of the constituent links by 
comparison with the axial maps. The Q = 6 field constructed using the rational map (13.71) 
with n = 2 is displayed in Figure [2b. As can be seen from this figure, the field is qualitatively 
similar to the minimal energy Q = 6 soliton, which is of the type £^2- More examples of 
rational maps associated to linked configurations will appear in the following section. 



4 Higher charge solitons 

In this section the results of a large number of energy minimization simulations are 
presented. Hopf charges up to sixteen are studied, with initial conditions consisting of a 
variety of links and knots, created using the rational map ansatz discussed above. 

Table 1 summarizes the simulations performed, by listing the type of each initial condition 
together with the type of the resulting final solution, obtained from the energy minimization 
algorithm. In Table 2 the types and energies of the known lowest energy solutions, plus 
some local energy minima, are presented for 1 < Q < 16, ordered by increasing energies. 
Also given in Table 2 are the energies divided by Ward's conjectured bound, E/Q 3 ^, from 
which it can be seen that the minimal energy solution is consistently around 20% above 
the conjectured bound. The values of E/Q 3 / 4 are plotted in Figure [3] using a notation that 
identifies the different torus knots as follows. White circles denote configurations which do 
not contain knots, that is, either unknots or links. The symbols denoting knots are triangles 
for /C 3i2 , diamonds for /C 5j2 , squares for /C 4)3 and stars for /C 5j3 . Configurations denoted by 
black circles will be discussed shortly. 
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Table 2: Solution types and energies. 
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Figure 3: The ratio of the energy to the conjectured bound, that is, E/Q 3 ^ 4 , as a function of 
Hopf charge Q for a variety of solutions: unknots and links (white circles); K,^ (triangles); 
£■5,2 (diamonds); JC^ (squares); K.^ (stars); links which are not resolved (black circles). 



The present study has produced 26 new solutions with 8 < Q < 16. The position curves 
for each of these is displayed in Figure HI but for clarity the linking curves are not shown, 
although they have been examined to confirm the correct linking number identifications. 
Each plot is labeled by its charge and type, with energies increasing first from left to right 
and then top to bottom. 

First of all, consider trefoil knots, that is, solutions which have the form As men- 
tioned earlier, the Q = 7 minimal energy soliton has this form and is obtained from the 
related rational map initial condition. This is encoded in Table 1 as the process ^3,2- 
The entry to the right of this one in Table 1 reveals that the same trefoil knot solution is 
also obtained from the linked initial condition C 2 's ■ 

There are no trefoil knots with Q < 7. Evidence supporting this is presented in Table 
1 for Q = 5 and Q = 6, where it is seen that in both cases initial conditions of the form 
/C3 2 result in the linked minimal energy solutions, which are C\;\ and £ 2 ' 2 respectively. A 
reasonable interpretation of these results is that for Q < 7 the number of twists T = Q — 3 
is too low, given the preferred length of the soliton in a trefoil knot arrangement. In other 
words, the twist per unit length is too low to be an energetically efficient distribution of the 
Hopf charge between crossing and twisting. A more detailed discussion of this aspect will 
be given later, when general torus knots will be considered. Note that, in particular, there 
is no trefoil knot solution in which there is no twist, which corresponds to Q = 3, and was 
the original suggestion for a knot soliton [5] . 
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Figure 4: The position curves for a variety of solutions with Hopf charge 8 < Q < 16. Each 
plot is labeled by its charge and type. 
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Table 1 shows that there are also trefoil knot solutions for Q = 8,9, 10, 11 and in each 
case they are obtained from rational map initial conditions of the type For Q = 9 the 
trefoil solution is also obtained from a variety of other initial conditions, including links and 
the torus knot JC^] this last process — > £3,2 provides an example of knot transmutation, 
where Solomon's seal knot deforms into a trefoil knot. In Figure H] a comparison of the plots 
8^-3,2; 9/C3 5 2, 10X^2, H^3,2) emphasizes that the conformation of the knot is important, with 
each of these trefoils having a very different structure to the others. In particular, as the 
number of twists increases the knot increasingly contorts and bends, in a manner similar to 
that seen in the axial solutions A Ut i, where for n = 3 the response to the twist is to lower 
the energy by bending to break the axial symmetry. As mentioned earlier, for n = 4, 5 there 
are similar bent solutions, which are local minima, where the deformation increases with n. 
Trefoil knots therefore appear to follow a similar pattern, with a minimal number of twists, 
that is 4, required for existence and increasingly deformed solutions existing with slightly 
larger twists than the minimal value. 

The triangles in Figure [31 and the associated energy values in Table 2, show that the 
ratio of the energy to the bound steadily increases for the trefoil solutions as the charge 
increases. It is therefore not surprising that for a large enough charge, which happens to be 
Q = 12, a trefoil solution fails to exist. Table 1 confirms that a Q = 12 initial condition of 
the type K^^ produces a solution which is not a trefoil knot, and will be discussed later. The 
triangles in Figure [3] also confirm that only for the minimal value Q = 7 is the trefoil knot 
the global minimum energy solution, and in all other trefoil is only a local energy 

minimum. However, for Q = 8 the trefoil knot energy is extremely close to the lowest energy 
found, which is a linked solution £33. In fact the energies differ by less than 0.1%, which 
is probably smaller than the accuracy of the computations, so in this case it is difficult to 
make a definitive statement regarding which (if any) is the global minimum. 

In Ref. [1] a linked solution of the type £22 was reported, but this appears to be an artefact 
of that study being at the limit of computational feasibility for computing resources available 
at that time. A configuration of this type does appear during a relaxation procedure, but 
after further relaxation it changes its type. Moreover, an initial condition of the type £2 2 
can be constructed using the rational map 

74 72 72 
W = 1 = ^ + ^ (4 1) 

and, as seen from Table 1, the relaxation yields a solution of a different type, namely the 
linked solution £3^3. It therefore now seems unlikely that a charge eight stable solution exists 
of the type £%\- Note that the rational map ( 14. ip describes a link in which both components 
are linked twice with each other, and this is because the irreducible factors of the denominator 
are terms like {Z\ — Z ), rather than terms like {Z\ — Z ), where the latter corresponds to 
only single links between any two components. The partial fraction decomposition in (14.11) 
reveals that each component of the link has charge two, by once again reading off the power 
of the numerator. 

For a given fixed charge there is nothing to prevent multiple trefoil knot solutions which 
differ in their conformations and energies, but no evidence for this phenomenon has been 
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found. For example, Table 1 shows that the Q = 9 trefoil knot is produced from (at least) 
four different types of initial condition, but in each case the resulting trefoil solution is the 
same one with an identical conformation. This is the case for all the solutions described in 
this paper, that is, when a particular solution is obtained from various initial conditions it is 
always in the same conformation, so its charge and type are sufficient to distinguish it from 
any other solution. 

As the trefoil knot solutions for Q = 9, 10, 11 are clearly only local energy minima then 
it remains to describe the candidates found for the global minima for these charges. In fact 
for all three of these charges, and also for Q = 12, the lowest energy solutions found have a 
similar form, which is a link with three components. The simplest case, Q = 9, is displayed 

2 2 2 

in Figured] as plot 9£i'i\, and as this label suggests there are three components to the link, 
each of which has charge one and links once with each of the other two components. An 
initial condition of the type £i'i'i with Q = 9 is obtained from the rational map 

Z 3 

W = — (4 2) 

and quickly relaxes to the minimal energy solution with the same type. The initial condition 
has a more planar arrangement of the three components than the final solution, and the 
initial cyclic C3 symmetry is also broken, as can be seen from plot 9£^'i, in Figure HI Note 
from Table 1 that this 3-component link is also obtained from a knotted initial condition of 
the type JC ij3 , but that out of the six Q = 9 initial conditions used, four lead to the higher 
energy trefoil knot solution. This is evidence that supports the fact that minimal energy 
solutions may not be the easiest to find; hence the need to employ the variety of starting 
configurations used in this study. 

The minimal energy solutions with Q = 10, 11, 12 are similar to the minimal Q = 9 solu- 
tion and correspond to increasing one, two and finally all three, of the charge one components 
to charge two. In other words they have the forms £1'? 2, £1 2^2 an d £2 2I2 respectively. These 
solutions are displayed in Figure H] and in each case there is an associated rational map of 
the same type. Table 1 shows that they can all be obtained from the relaxation of certain 
torus knot initial conditions. 

So far the only knot solutions discussed have been trefoil knots. As the torus knot 
with the lowest crossing number after the trefoil is the (5, 2)-torus knot, then this is the 
most likely candidate to appear next as a soliton solution. Table 1 reveals that an initial 
condition of the type £5,2 does not yield a solution of this type for Q = 8, 9, 10, but it does 
for Q = 11. The Q = 11 solution is displayed in Figure 0] as plot HK,^- It is only a local 
energy minimum, despite the fact that it is lower in energy than the Q = 11 trefoil knot. 
There are also solutions of the type for Q = 12, 13, 14 and the pattern mirrors that of 
the trefoil knots, in that the ratio of the energy to the bound steadily increases with the 
charge; see the diamonds in Figure [3j 

At Q = 12 a solution of the type /C 4j3 appears (see Figure 0] plot 12/C4 ;3 ) and its energy 
is just above that of the minimal energy £^'2 2 link, though substantially lower than that of 
the /C 5)2 solution. Solutions of the type /C 4i3 also exist for Q = 13 and Q = 14, and for these 
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two charges they are the lowest energy solutions found. The excess energy above the bound 
is also reasonably low, so these two solutions are good candidates for the global minima at 
Q = 13 and Q = 14. The squares in Figure [3] denote the ratio of the energy to the bound 
for K,^ solutions, from which it can be seen that (unlike the other torus knots) the lowest 
charge at which this knot appears is not the one that is closest to the conjectured bound. 

At Q = 13 a solution exists (see Figure 0] plot 13Xi 3 ) which is denoted by X 13 because 
there is no unambiguous interpretation as a particular link or knot. This is due to the fact 
that a definition of a link or knot requires the position curve not to self- intersect, but there is 
nothing to prevent this in the field theory. For the X i3 solution then either the position curve 
self-intersects exactly, or there are parts of the curve that are so close together that they can 
not be resolved with the numerical accuracy currently employed. Some of the lower charge 
solutions presented earlier may also appear to self-intersect, for example the Q = 5 solution 
displayed in Figure [TJ but in all these lower charge solutions the apparent self-intersection 
can be resolved by reducing the thickness of the tube plotted around the position curve, 
together with a careful consideration of continuity. 

If there are no self-intersection points in the solution Xi 3 then there are two possibilities 
for the configuration type, depending on how this self-intersection point is resolved. The 
first possibility is the link C^ 2 i an d the second is a link with two components, where one of 
the components is a trefoil knot with five twists and the second component is approximately 
axial with a single twist. The latter resolution is essentially an intertwining of the charge 
eight trefoil knot solution with the charge one solution. 

The lowest energy solitons found for Q — 15 and Q — 16 are denoted by X\§ and 
Xiq respectively, and also have points which can not be distinguished from self-intersection 
points, as seen from plots 15Ai5 and I6A46 in Figure HI In each case there is a resolution into 
a link with two components where one of the components is a trefoil knot. The conformation 
of both the knot and unknot components in these cases strongly suggests that this is the 
correct resolution, if indeed one is required. The Q — 16 solution is the easiest to identify 
and consists of the charge eight trefoil knot intertwined with the minimal energy charge 
two solution, where the conformation of both components is very similar to that of the 
components in isolation. A field with the correct qualitative behaviour to describe this 
resolution is given by the rational map 



Initial conditions generated using this rational map quickly leads to the solution X\q under 
energy relaxation. 

The solution A45 is similar to the solution Xi$, except that the resolution involves the 
charge seven trefoil knot, rather than the charge eight trefoil knot. The energies of the 
solutions of type X n are represented by black circles in Figure [3j 

Note that Q = 16 is the lowest charge in which a link with four components each linking 
all the others might exist. However, as seen in Table 1, an initial condition of the type £i'x'i'i 
relaxes to the X w solution. 



W 



Zf + Z^Zq + Z\Zq 
ZoZf + Zq 




+ 



(4.3) 
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Figure 5: String length (circles) as a function of charge Q for several solitons, together with 
the expected growth behaviour ^Q 3 ^ A (curve). 



The final torus knot solitons presented in this paper are of the type £5,3 and exist for 
Q = 14 and Q = 15, but in both cases these solutions are only local energy minima. The 
energies of these solutions are represented by stars in Figure [31 

In summary, it has been shown that there are a variety of torus knot solitons at various 
Hopf charges. For example, at Q — 14 three different torus knot solutions have been obtained. 
Most of the knot solitons found are only local energy minima, but in some cases they are 
good candidates for the global minimum. 

For a string-like soliton it is expected that one of the contributions to the energy should 
have an interpretation in terms of the string length. As the soliton energy grows like Q 3 / A 
then this suggests that the string length should have a similar growth. In Figure the string 
length is plotted (circles) for several solitons, including examples of both knots and links. 
Also shown is a curve with the expected power growth, 7Q 3 / 4 , where 7 has been computed 
by a least squares fit to be 7 = 7.86. This plot demonstrates a reasonable agreement with the 
expected Q 3 ^ dependence of the string length. The greatest discrepancy from the expected 
string length occurs for the Q = 15 link C^W, which is much longer than expected. However, 
note that this solution has an abnormally large ratio of linking to twist number, with the 
contribution to the Hopf charge from linking being four times that from twisting. The 
solution with the next largest ratio of linking to twist is the Q = 9 link C 1 ' 1 \, where linking 
contributes twice as much as twisting, and the length of this solution is also a bit above the 
expected value. It therefore appears that solutions which link much more than they twist 
require an increased string length to minimize energy. 
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Recall from the earlier discussion of knot solitons, that there appears to be a critical twist, 
below which knot solitons of a particular type do not exist. A twist below the critical value 
seems to be an energetically inefficient distribution of the Hopf charge between twisting and 
crossing. For twists just above the critical value a knot solution continues to exist but as the 
twist increases further then again the distribution between twisting and crossing becomes 
inefficient, but this time due to too much twisting. 

A naive approximation is to assume that for all solutions there is a universal optimum 
twist per unit length. Taking this optimum value from that of the Q — 1 soliton, and using 
the fact that the string length grows like Q 3 ^ 4 , this suggests that an optimum value for the 
number of twists is approximately T ps Q 3 ^ 4 . For a knot with crossing number C this gives 
C = Q — T Q — Q 3 / 4 . Using this assumption, a knot with crossing number C is predicted 
to be energetically efficient at an integer charge near to the real value Q*, which is defined 
as the solution of the equation — Q 3 J 4 = C. For example, for the trefoil knot C = 3 and 
therefore Q* = 7.56, which roughly agrees with the Q = 7 and Q = 8 trefoils being closest to 
the energy bound. For C = 5 then = 11.07 and this explains why the Q — 11 knot K,^ 
is closest to the bound. For C = 8 then Q* = 16, so the naive approximation overestimates 
the most efficient value, which is seen to be Q = 13 for the /C4 3 knot. There are clearly more 
subtle effects at work than the simple naive constant twist rate assumed in this calculation, 
but it does seem to produce numbers which are in the right ballpark, suggesting that it 
captures some qualitative aspects of knot energetics. 

5 Conclusion 

The results presented in this paper reveal that there are many low-energy knot solitons 
of various types, for a range of Hopf charges, together with increasingly complicated linked 
solitons. The qualitative features of these solutions can be replicated using an ansatz involv- 
ing rational maps from the three-sphere to the complex projective line, and futhermore this 
provides a good supply of initial conditions for numerical relaxation simulations. The com- 
plicated nature of the problem means that it is difficult to be certain that the global minimal 
energy solitons have been found, but certainly some good candidates have been presented, 
whose energies agree well with the expected values based on an earlier conjectured bound. 

The rational map ansatz allows the construction of any torus knot initial condition, but 
it is not clear how to extend this to non-torus knots. All the knot solutions found so far 
are torus knots, and it is unknown whether the lack of non-torus knot solitons is a result 
of not having suitable initial conditions, or whether there is some energetic reason to favour 
torus knots. For example, the figure eight knot has four crossings and hence a soliton with 
this form might be expected with a charge around Q — 9. As seen from Table 1, quite a few 
initial conditions have been employed for Q = 9 and Q = 10, and no such solution has been 
found. It remains an open problem to understand the absence (or otherwise) of non-torus 
knots. 

Finally, given that there are many knot and link solitons then it would be useful if there 
was an approximate string model that could, even qualitatively, reproduce the field theory 
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results. It has been shown that the string length has the expected behaviour with Hopf 
charge, so it is plausible that there may be an approximate description based on a string 
energy that includes contributions from properties of the string such as its length, twist and 
writhe, together with relevant interaction terms. 
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